Parallel magnetic field tuning of valley splitting in AlAs two-dimensional electrons 
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We demonstrate that, in a quasi- two-dimensional electron system confined to an AlAs quantum 
well and occupying two conduction-band minima (valleys), a parallel magnetic field can couple to 
the electrons' orbital motion and tune the energies of the two valleys by different amounts. The 
measured density imbalance between the two valleys, which is a measure of the valley susceptibility 
with respect to parallel magnetic field, is enhanced compared to the predictions of non-interacting 
calculations, reflecting the role of electron-electron interaction. 

PACS numbers: 



The physics governing the valley splitting in multi- 
valley two-dimensional electron systems (2DESs), such as 
in Si field-effect-transistors, has been of interest for some 
time [3, 01 ■ It is attracting renewed attention in 2DESs 
in Si, AlAs, and graphene [U, HI, both as a fundamental 
problem, and also because of the possibility that manip- 
ulating the electron valley degree of freedom might lead 
to novel (" valleytronics" ) devices [!, E| 0] ■ Traditionally, 
the valley energies have been controlled via strain, con- 
finement and electric field 0, 0, H 0| ■ In this Letter we 
demonstrate how a magnetic field (-B||) applied parallel 
to the 2D plane can also break the valley degeneracy 
and shift the valley energies of a 2DES with finite layer 
thickness. In such a system, B|| couples to the electron 
orbital motion and deforms the electron wave function 
in the confinement direction. To first order, this defor- 
mation increases the valley energies (diamagnetic shift) 
and, to second order, it increases the effective mass in 
the direction perpendicular to B» [ljj, El, 12, EH, Ell- 
in our AlAs samples, where two valleys with anisotropic 
Fermi contours are occupied, when B\\ is applied along 
the major axis of one the valleys and perpendicular to 
the other valleys' major axis, it shifts the valleys' en- 
ergies by different amounts. Remarkably, the measured 
energy shift and the resulting valley density imbalance 
are much larger than simple calculations would predict, 
signaling a clear enhancement of the -By-induced valley 
spitting, likely due to electron-electron interaction. (The 
parameter r s , defined as the average intcr-electron spac- 
ing measured in units of the effective Bohr radius, ranges 
from 6.4 to 9.8 for our samples.) 

Figure 1 summarizes our experimental setup and mea- 
surements. We studied four samples; here we focus on 
two samples where the 2DES is confined to either an 
11 nm-widc or a 15 nm-wide AlAs quantum well (QW), 
grown using molecular beam epitaxy on a semi-insulating 
GaAs (001) substrate. The AlAs well is flanked by Al- 
GaAs barriers and is modulation-doped with Si [l5j. In 
these samples the electrons occupy two conduction-band 
valleys with elliptical Fermi contours as shown in Fig. 
1(b) [3], each centered at an X point of the Brillouin 
zone, and with an anisotropic effective mass (longitudi- 
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FIG. 1: (Color online) Schematic diagrams showing: (a) Ex- 
perimental setup; (b) and (c) Valley occupations at zero (light 
green) and finite (dark green) parallel magnetic fields; (d) and 
(f) Fermi contour distortions due to parallel magnetic field 
for the [100] and [010] valleys; (e) and (g) Density-of-states 
(DOS) for the [100] and [010] valleys before (light green) and 
after (dark green) the application of parallel magnetic field. 



nal mass m; = 1.05 and transverse mass m t = 0.20, in 
units of free electron mass, m e ) We denote these two 
valleys according to the direction of their major axis: 
[100] and [010] (Fig. 1(b)). In our experiments B\\ is 
applied along [100] as schematically shown in Fig. 1(a). 
Since the two valleys have different orientations with re- 
spect to -By, the diamagnetic shift and the effective mass 
enhancement are different for the two valleys, as illus- 
trated in Figs, l(d-g), causing an electron transfer from 
the [100] to the [010] valley (see Fig. 1(c)). This density 
imbalance (An) caused by B» can be countered by ap- 
plying symmetry-breaking strain e = e[ 10 o] — £[oio], where 
£[ioo] and ejoio] are the strain values along the [100] and 
[010] directions. In our study, we monitored the sample 
resistance vs. e at different values of -By to determine 
An as a function of Bu . In order to apply tunable strain 
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our samples were glued to apiezoelectric stack actuator 
as shown in Fig. 1(a) 0, 17 1. The measurements were 
performed in a 3 He cryostat with a base temperature of 
0.3 K. The system was equipped with a tilting stage, al- 
lowing the angle 9 between the sample normal and the 
magnetic field to be varied in situ. 

Before presenting the experimental data, we will first 
outline the theoretical formalism that describes this val- 
ley splitting. Using the same approach as in Refs. 
and hll, we start with a simple, single-particle effective 
mass Hamiltonian in three-dimensions: 
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where p Xl p y , p z are momentum operators, m 
are effective masses in x, y, z directions, q is the electron 
charge and V(z) is the confinement potential in the z 
direction which is assumed to be a QW with 210 meV 
barrier heights and 11 nm or 15 nm of well width for 
our two sample structures [l(|. We use the gauge A = 
(0, -B||2i, 0). Solutions to the above Hamiltonian are plane 
waves in x and y directions. Substituting the plane wave 
solutions leads to the one-dimensional Hamiltonian: 
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where the last two terms containing S11 can be described 
as perturbations to the Hamiltonian H = p 2 z /2m z + 
V{z). To first order, the ground state energy is shifted to 
higher values because of the second perturbation term, 
(qB\\z) 2 /2m y . To second order, the effective mass in the 
y direction (h 2 j \d 2 E j dk 2 )) becomes enhanced since the 
first perturbation term (hk y qB^z/m y ) is linear with k y . 
Eventually, the total energy for the ground state can be 
written as follows: 
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where E is the ground state energy of H , AE is the 
ground state energy shift, and m y is the enhanced effec- 
tive mass in the y direction. 

Initially, the two valleys have the same solutions to 
Hq because of their same mass in the z direction. The 
crucial observation is that the perturbation terms in Eq. 
(2) depend on m y which is different for the two valleys. 
The [100] valley has a smaller mass in the y direction and 
therefore feels a stronger perturbation at a given B^ . Its 
energy is shifted more compared to the [010] valley, caus- 
ing a splitting with B|| . The quantity we measure exper- 
imentally is An caused by B», and contains two terms: 
one, the difference in the ground state energies and the 
other, the difference in the density-of-states caused by 
the mass enhancement. To second order in perturbation, 
An can be written as: 
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FIG. 2: (Color online) (a) Resistance vs. strain traces for 
an 11 nm-wide AlAs QW taken at v = 3 for different par- 
allel magnetic fields. The "balanced point" for a given B\\ 
is marked by a closed circle on the corresponding trace, (b) 
Landau level fan diagram for Bu = 3.4 T as a function of e. 



where a = {l/2)q 2 B 2 ((z 2 } 



(z) 2 ) and Ei is the first ex- 



cited state energy of Hq- Although second order pertur- 
bation theory gives analytic answers, in our calculations 
we numerically solved the Schrodinger's equation using 
a finite difference method; i.e., the values we report are 
numerically exact solutions to Hid in Eq. (2). 

Experimentally, we determine An vs. B\\ by moni- 
toring the sample resistance (R) as a function of e in 
tilted magnetic fields. Examples of such piezo-resistance 
traces are shown in Fig. 2(a). Each trace was taken at a 
fixed 9 and magnetic field so that the 2DES remains at 
a fixed Landau level (LL) filling factor (^=3 in the case 
of Fig. 2(a) data). With applied e, the LLs for the [100] 
and [010] valleys cross each other, as the fan diagram in 
Fig. 2(b) indicates. R exhibits minima as the Fermi en- 
ergy (Ep) passes through consecutive energy gaps, and 
maxima as it coincides with the LL crossings. The fan 
diagram in Fig. 2(b) is drawn for a fully spin polarized 
system; this is indeed the case for the traces shown in 
Fig. 2(a) because of the large Zeeman splitting due to 
the high magnetic field. From the piezo-resistance traces, 
the "balanced point," i.e., the strain at which the valleys 
are equally occupied at a given B\\ can be measured by 
following the resistance peaks indicated by closed circles. 
(The balanced point at B\\ — 0, which corresponds to 
e = 0, is determined experimentally from R vs. e sweeps 
at even fillings as has been detailed in Ref. 17J-) The 
data of Fig. 2 clearly show that the valley energies are 
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split with the application of B\\ . 

In Fig. 3(a) we show a set of R vs. e traces taken at 
a fixed 9 for v = 4, 5 and 6. Again, in all these traces, 
the 2DES is fully spin polarized because of the very large 
Bn. The traces are periodic in e, implying that they 
are consistent with a simple, linear LL fan diagram, an 
example of which is shown in Fig. 6(b) for v = 6. By 
associating the R maxima with the LL coincidences, we 
can determine the valley polarization, Py , defined as the 
difference between the [010] and [100] valley populations 
divided by the total 2DES density (i.e., An/n), at each 
coincidence. We therefore obtain a direct measure of Py 
vs. e which we plot, for v = 6, as black squares in Fig. 
3(c); note that Py is equal to -1/6, 1/6, 3/6 and 5/6 for 
v = 6 coincidences from left to right. The lines in Fig. 
3(c) provide two pieces of useful information. First, their 
intercepts with the Py axis provide direct measures of 
Py, or cquivalently, An, for given values of Bp In the 
inset to Fig. 3(c), we show a plot of An vs. Bp Note that 
the plot is approximately linear, consistent with what we 
expect from Eq. (4). Second, the slopes of the lines in 
Fig. 3(c) give a measure of the valley susceptibility with 
respect to strain [T3|, i.e., the rate of change of An with 
e; note that the slopes are independent of v for a fixed 
density, consistent with the measurements of (l7l | which 
were done in the absence of parallel field. 

To summarize our results, in Figs. 4(a) and (b) we 
plot the measured An vs. e and vs. Bp respectively, 
for the 11 nm-wide AlAs QW at three different densities. 
The data in Fig. 4(a) were obtained from lines such as 
those shown in Fig. 3 (c), but after subtracting their in- 
tercepts with the An axis. Figure 4(a) data essentially 
represent An vs. e in the absence of Bn\ this can be 
verified from comparison of the data with measurements 
reported in Ref. [L7| which were done at 9 = 0. In Fig. 
4(a) we also include a plot of An vs. e expected based on 
the band parameters, i.e., ignoring electron-electron in- 
teraction and using the simple relation An = eE2m/2'KTi 2 
where m = y/mjmt =0.46 and £^=5.8 eV is the defor- 
mation potential for the AlAs X-point conduction band 
minimum [l6| . In Fig. 4(a), it is clear that the response 
of the system to strain is two to three times enhanced 
compared to what the non-interacting (band) parameters 
predict, and that the enhancement is stronger at lower 
densities (larger r s ). This observation confirms the en- 
hancement of valley susceptibility with respect to strain 
(Xe), defined as the rate of change of An with e, originally 
reported in Ref. [l?] ]: it is similar to the enhancement of 
the spin susceptibility observed in similar samples and re- 
flects the role of electron-electron interaction jl7 , HI, 3] • 



Our measured An caused by B^ are plotted in Fig. 
4(b) vs. Bjj for the same three densities as in Fig. 4(a). 
Again, for comparison, we also show the results of the cal- 
culations based on the non- interacting picture, i.e., using 
Eq. (4). As can be seen in Eq. (4), there is a slight de- 
pendence of An on n and this is why there arc three lines 
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FIG. 3: (Color online) (a) Resistance vs. strain traces for the 
same sample as in Fig. 2 but at a higher density, taken at 
different fillings for a fixed tilt angle, 9 — 74.5°. (b) Landau 
level fan diagram as a function of e for v — 6. (c) Valley 
polarization (Pv) as a function of e at different fillings. The 
intercepts of the lines with the y axis give Pv (or alternatively, 
An) for the corresponding values of B\i . Inset shows An as a 
function of B» . 



in Fig. 4(b) representing the non-interacting calculations 
pol ]. Similar to the strain case of Fig. 4(a), the sys- 
tem's response to is about two to three times stronger 
than the non- interacting calculations predict. This is a 
noteworthy result as it demonstrates that the interacting 
2DES responds to two very different stimuli (strain and 
Bu) in a very similar fashion. Defining a new quantity, 
XB||, as the valley susceptibility with respect to Bp we 
see from the inset to Fig. 4(b) that Xe an d Xb« have very 
similar magnitudes [2l| . 



In Fig. 4(c) we show the summary of -By-induced An 
for a second sample, a 2DES confined to a 15 nm-wide 
AlAs QW. The data overall arc quite similar to the 11 
nm-wide QW data. The larger well- width of the 15 nm- 
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FIG. 4: (Color online) (a) Summary of density imbalance as function of e for three densities for the 11 nm-wide QW. (b) and 
(c) Summary of density imbalance as function of By for three densities for the 11 nm-wide and 15 nm-wide AlAs QWs. Insets 
show the valley susceptibilities, normalized to their band values, with respect to two variables, e and Sjf. 



wide sample predicts a larger An for this sample (at a 
given density) because of a larger spread of the wave func- 
tion in the z direction (see Eq. (2)). The deduced, nor- 
malized susceptibilities plotted in Fig. 4(c) inset, how- 
ever, appear to be slightly smaller than those for the 11 
nm-wide sample. It is likely that this is also related to the 
larger electron layer thickness: Although there have been 
no systematic and detailed studies of the dependence of 
Xe enhancement on electron layer thickness, it has been 
reported that the spin susceptibility of a thicker electron 
system is less enhanced compared to a thinner system 
with otherwise the same parameters Til 22|. 



Our results demonstrate that a magnetic field applied 
parallel to an AlAs 2DES with finite layer thickness can 
split the energies of the two in-plane valleys. The split- 
ting originates from the coupling of the parallel field to 
the orbital motion of the electrons. Our measurements of 
the density imbalance due to strain and parallel magnetic 
field show similar enhancements of the splitting com- 
pared to calculations which are based on non-interacting 
electrons with band parameters. These enhancements 
are more pronounced at lower densities (larger r s ). The 
results suggest that in an interacting 2DES, the valley 
splitting is enhanced by interaction, independent of the 
physical parameter that induces the splitting. 
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